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Abstract 

We investigate the critical dynamics of 0(iV)-symmetric scalar field theories to determine the 

6-i 

D . critical exponents of transport coefficients as a second-order phase transition is approached from 



the symmetric phase. A set of stochastic equations of motion for the slow modes is formulated, 
and the long wavelength dynamics is examined for an arbitrary number of field components, N, 
in the framework of the dynamical renormalization group within the e expansion. We find that 
for a single component scalar field theory, N = 1, the system reduces to the model C of critical 
dynamics, whereas for N > 1 the model G is effectively restored owing to dominance of O(N)- 
symmetric charge fluctuations. In both cases, the shear viscosity remains finite in the critical 
region. On the other hand, we find that the bulk viscosity diverges as the correlation length 
squared, for N = 1, while it remains finite for N > 1. 
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I. INTRODUCTION 



In recent decades transport coefficients of the quantum chromodynamics (QCD) have 
attracted much interest in the context of Relativistic Heavy Ion Collider (RHIC) experiment, 
which aims at creating and studying a quark-gluon plasma. One of the interesting findings 
emerging from the experimental program at RHIC, the large elliptic flow v 2 observed in high 
energy non-central collisions, implies that the spatial anisotropy of the initial state created 
in the collision is efficiently converted during the expansion to a transverse momentum 
anisotropy of the observed hadrons These experimental results are well described by 

ideal hydrodynamics with vanishing viscosity {4-6]. Thus, the large elliptic flow observed 
in such collisions implies that the matter created in collisions behaves as an almost perfect 
fluid. 

Although the transport coefficients in viscous hydrodynamics are phenomenological pa- 
rameters, they can, in principle, be computed from a microscopic theory. Since the shear 
viscosity, one of the transport coefficients, has a direct influence on the elliptic flow, the 
experimental results have triggered numerous theoretical efforts to unravel its behavior as 
a function of thermodynamic variables. In general, these are performed in the framework 



of 



dnetic theory, e.g. using the Boltzmann equation, applied to effective theories of QCD 



14j and to perturbative QCD 15l-ll9|. Furthermore, some results on the temperature 



dependence of the transport coefficients have been obtained in lattice simulations |2Q|-|26|. 

The results of the RHIC experiments have motivated recent work on a field theoretical 
approach to evaluate transport coefficients. The O(N) scalar field theory offers a testing 
ground for developing computational methods before facing the complications of a full QCD 
calculation. The scalar field theory has in fact long been studied as a prototype theory in 
many contexts of physics. Moreover, for N = 4, the O(N) model serve as a low energy chiral 



271 ] . A general Lagrangian density for the O(N) scalar 



effective theory for two flavor QCD 
field theory is given by 

C = \m?-\rti-\u{<Pl)\ (1) 

where r is the mass parameter, u is the coupling constant, and the implicit summation over 
% runs from % — 1 to N. Recently, the dynamical properties of the O(N) scalar theory, in 
particular behavior of the transport coefficients, formulated microscopically in the Green- 
Kubo-Nakano linear response theory, have been explored in several theoretical studies. 



The shear viscosity, 77, of the scalar field theory was first studied by Hosoya et. al. 28( and 



Jeon and Yaffe 



29 



by Aarts and Resco 



30 
31 



in thermal field theory. Later on, the large N behavior was examined 
32| . These calculations demonstrated that rj is an increasing function 



of temperature, T. At high temperatures, 



V 



u 2 



(2) 



The cubic power in temperature can be understood on dimensional grounds, and the factor 
N 2 is attributed to the scaling of the coupling constant u with 1/N. The inverse power of the 
coupling constant in Eq. ([2]) implies that the shear viscosity is a non-perturbative quantity. 
The precise numerical factor in Eq. (T5|) can be obtained by a resummation of ladder type 
diagrams. In Refs. [29 



301 ], it was found that the ladder resummation is equivalent to the 
linearized Boltzmann equation with a thermal mass term. Some systematic approaches for 
computing higher order corrections are presented in Refs. 
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the transport equation are discussed in Refs. [36L l37|. 



351 ] . and relevant issues on 



In the present work, we discuss the critical behavior of the shear viscosity and other 
transport coefficients in the O(N) scalar field theory. As demonstrated by Wilson using 
the renormalization group approach, there is a second-order phase transition in the O(N) 
d theory. Kinetic approaches employed for computing the transport coefficients 



scalar fie 
(see e.g. 



38 



and a discussion in Ref. 



40] ) rely heavily on Boltzmann-like approxima- 



tions, which take only the single particle distribution into account and neglect higher order 
correlations. Although these correlations, may be unimportant far from the critical point, 
they play an important role in the critical region. 

In our study of the critical transport properties, we employ the dynamical renormalization 
group (DRG) combined with the epsilon expansion * developed by Hohenberg and Halperin 
(for a review, see Ref. j^). Within this approach we examine the scale evolution of a 
stochastic equation of motion, which describes the critical dynamics of slow modes. These 
include fluctuations of the order parameter and of conserved quantities, which are relevant 
variables when addressing the long- wave length behavior of the system near the critical point. 



* An alternative non-perturbative approach to the epsilon expansion would be a direct application of the 
functional renormalization group (see e.g. Ref. 4l[ for a review) to a quantum-field model constructed 



to b e eq uivalent to the stochastic equations of motion 
Ref. (44 1. 
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431. This method was tested for model A in 
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Since the transport coefficients are obtained from the corresponding response functions by 
taking the limit of both frequency and momentum to zero, they characterize the dynamics 
of the system in the low energy limit. 

In analogy to the static case, the flow equations for transport coefficients derived from 
the DRG admit non-trivial fixed points, from which the dynamical critical exponent, z, and 
the dynamical scaling relations can be derived. The dynamical critical exponent, z, defines 
the characteristic frequency of the most relevant slow mode cu ~ k z , and the scaling relations 
link the singular contribution to the transport coefficients. From these properties one can 
deduce the singular behavior of the transport coefficients, in particular whether they diverge 
or remain finite at the critical point. Based on the universal behavior, i.e. on the dynamical 
critical exponents and scaling laws, one identifies each system with a dynamical universality 
class. In contrast to the static case, the dynamical universality class is governed not only 
by the dimensionality, locality, and the symmetries of the system under consideration, but, 
in addition, by the properties of the relevant slow modes. Thus, the conservation or non- 
conservation of an order parameter, and the existence of mode-mode couplings among the 
slow modes affect the dynamical universality class. Therefore, even if two systems belong to 
the same static universality class, their dynamic universality class ' may be different *. In 



t In what follows, we will frequently refer to the universality classes that were defined in Ref. 45 1. Here we 
provide some properties of the relevant universality classes: 



Model 


Slow mode(s) 


Dynamical critical exponent in d=3 


A 


Non-conserved field 


z A = 2 + 0.7621?7' 


B 


Conserved field 


zb = 4 - rf 


C 


N-component non-conserved field 


z c = 2 + a/v for TV = 1 




coupled to one component conserved field 


z c = 2 + 0.7621??' for N > 1 


H 


Conserved field coupled to 
conserved transverse vector field 


z H = 4- 18/19 


G 


N-component non-conserved field 
coupled to N(N — l)/2-component conserved field 


z G = 3/2 
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Here a and v are the static critical exponents and rf is the anomalous dimension. 
* An example of such a situation is given by the models A and B of critical dynamics (see Refs. 
for further details). In the static case, both models, exhibiting Z(2) symmetry and belong to the static 
universality class of the Ising model in three dimensions. The dynamical universality class is, however, 
different. This difference, arises from the non-conservation (conservation) of the order parameter in 
model A (B), and results in different long wavelength behavior characterized by the dynamical critical 
exponent: za = 2 + const ■ rf and zb = 4 — rf where rf is the anomalous dimension. Another nontrivial 
example is the O(N) model for N = 1 (non-conserved order parameter) and model H (conserved order 
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this article, we determine the dynamical universality class of the 0(N) scalar field theory 
and show how the dynamical universality class depends on the number of components, N, 
and on the dimensionality, d. 

The paper is organized as follows: in the next section we identify the slow modes in 
the O(N) scalar theory, and construct an effective Hamiltonian for them. In section III we 
review the static universality classification of the theory, and show that a non-trivial fixed 
point exists. In section IV we introduce the stochastic equation of motion, which describes 
the dynamics of the slow modes in the critical region. We then implement the DRG to 
find the fixed points of the stochastic equations of motion , and determine the dynamical 
universality class. We close this section with a brief discussion of the critical behavior of 
the bulk viscosity. Section V is devoted to summary, discussion, and outlook. Details on 
the derivation of the stochastic equation of motion and on the calculation of the response 
function are given in two appendices. 

II. CONSTRUCTION OF THE EFFECTIVE HAMILTONIAN 

Before considering the dynamics of the theory, we have to build an effective Hamilto- 
nian for the slow modes A\_. The probability distribution for the modes Ai is given by the 
exponential of the Hamiltonian §, e~ n ^ Al ^. The effective Hamiltonian defines the static 
critical behavior of the theory, and will later on be incorporated in the equations of motion, 
from which we finally find the dynamical properties of the system close to the critical point. 
Although, the effective Hamiltonian, and the equations of motion for the slow modes have 
a microscopic origin, it is in general a very challenging problem to derive them starting 
from the microscopic Lagrangian. Therefore, in the present work, we formulate the effective 
Hamiltonian and the equations of motion on a phenomenological basis. The guiding prin- 
ciples in such a formulation are similar to those of Ginzburg-Landau theory. Note that, in 
our case, the slow variables in the Ginzburg-Landau Hamiltonian are all fluctuations, i.e., 

parameter) of critical dynamics. These models also share the same static universality class, while the 
dynamic universality classes differ. This implies a completely different behavior of quantities such as the 
shear viscosity, close to the critical point. In model H, the shear viscosity diverges, while, as will be shown 
in this article, it is always finite in the O(N) scalar field theory. 
§ Here, the prefactor in the exponent (fc^T) -1 is absorbed in the definition of the reduced effective Hamil- 
tonian. 
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deviations of variables from their equilibrium values. 

Candidates for the slow mode of the theory are the fluctuations of the order parameter 
<pi, the energy-momentum density, E and J, and the O(N) charge density, Qij. Owing to 
the symmetry = —Qji, there are N(N — l)/2 charges associated with generators of 
the O(N) group. The order parameter of the theory is not conserved, while the remaining 
variables (energy, momentum and O(N) charge) are conserved quantities. 

In the present work, we consider a system approaching the critical point from the sym- 
metric phase. In this case, it is straightforward to construct the effective Hamiltonian for 
the slow modes 



H = I d d x 



l r _i „, 1 1 



u+ + lotfE + -c, 1 e 2 + -J 2 + -xtQ% + n t 



(3) 



H+ = l{^) 2 + r jti + ^(ti)\ (4) 
% s = -fchi - J ■ H + {3E - HijQij, (5) 

where 1-L S is the source term, which is introduced for later convenience. We follow the con- 
vention that repeated indices imply summations, e.g., Q\j = \ X)ij=i QijQij- The effective 
Hamiltonian ([3]) includes all possible candidates for slow modes in an O(N) scalar field 
theory. 

Since the Hamiltonian includes up to quadratic terms in E, J and Qij, the original 
Hamiltonian density for the order parameter fluctuation, "H^, is recovered after integrating 
out these variables and performing a suitable redefinition of the couplings. This implies that 
the critical statics of the Hamiltonian % is the same as that of 

The coefficients of the Hamiltonian ((3]), are given by the static susceptibilities of the 
slow modes. Since the susceptibility of the momentum current J always remains finite, 
we have absorbed the coefficient of J 2 by a redefinition of the field J. The O(N) charge 
susceptibility, xo, also remains finite for zero net charge (i.e. zero chemical potential). In the 
case of Bose-Einstein condensation with a finite O(N) charge, however, xo diverges at the 
critical point. We do not consider this situation, but keep Xo explicitly in the Hamiltonian 
for later convenience. 

There are two contributions to the fluctuations of the energy density, 5E = T5S + 
hiStfii, where S is the entropy density^. Consequently, the static correlation with the order 



^ In the remainder of this section, we explicitly denote the fluctuation of a variable X by 5X , in order to 
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parameter fluctuation is given by 

(8E8<f>i) = T(5S5(j) j } + hjiSfaSfa). (6) 

In an O(N) symmetric system (no explicit symmetry breaking), hj = at the physical point. 
Therefore, the last term in Eq. ([6]) does not contribute to (5E5(j)i). The correlation (5S5<j)i) 
is nonzero at temperatures below T c , where the symmetry is spontaneously broken in a 
specific direction of the field. Thus, for this component of the field 0«: (SSSfa) ~ —<9t (</>«) ^ 
0. Moreover, this quantity diverges close to the critical point in the broken phase, since 
(5S5(pi) ~ t^ -1 , where t = (T — T c )/T c is the reduced temperature, and the critical exponent 
< 1/2. However, at temperatures above T c , the correlation function (SSSfa) vanishes due 
to symmetry in the absence of the external field hi. Indeed, since for t > and hi — > the 
order parameter scales like ~ t _7 /ij, which implies that (5S5<f)i) ~ t~ 1 ~ 1 h i Thus, in the 
effective Hamiltonian, there is no bilinear contribution of the form ~ <piE for t > 0. 
Now consider the autocorrelation function of the energy fluctuations 

(5E5E) = T 2 d T E, (7) 

which is proportional to the specific heat, C. Near the critical point, the singular part of 
C scales as ~ t~ a ~ ^ a l v where £ is the correlation length. The specific heat is related to 
the static susceptibility of the energy, C oc Xe(M = 0;T), up to some dimensionful factor. 
The sign and numerical value of the critical exponent a depends on the number of field 
components, N, and the dimensionality, d (see e.g. \^\)- 

In the effective Hamiltonian, we have dropped the spatial derivative terms, i.e. terms 
of the form (V m A;) 2 , for all fields A\ which turn out to be irrelevant for long wavelength 
physics, except for the order parameter. Consider for instance the term involving derivatives 
of the energy density, i.e. {^SJe\ . For negative a, the specific heat remains finite at the 
critical point. Hence, the coefficient of the E 2 term in the Hamiltonian scales as C~ l ~ £°. 
Using standard renormalization group arguments, one then finds that the derivative term 
{y E^j is irrelevant. Also for positive a, when the specific heat diverges as the term is 
irrelevant as long as ajv < 2. This inequality is in general satisfied, since ajv is small, 0(e), 
where e — 4 — d. In the case of interest, where the critical point is approached from the 

avoid ambiguities. 
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symmetric phase, i.e. T — > (T c ) + for /i = 0, all derivative terms of the conserved quantities 
are, by the same reasoning, negligible. Consequently, for static properties, contributions at 
the scale ~ A are due only to loop corrections involving fluctuations of the order parameter. 
The corresponding derivative term is relevant, yielding nontrivial contributions to the critical 
exponents through the nonzero anomalous dimension. 



III. CRITICAL STATICS 



A. Critical exponents and scaling hypothesis 
In tl 
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lis section we review the critical statics at continuous/second-order phase transitions 
48[. A general effective theory for the order parameter of a continuous phase transi- 
tion was developed by Landau. This theory provides a mean-field description of the phase 
transition. The Ginzburg criterion defines the region of applicability of the mean-field ap- 
proximation. Close to the critical point, in the critical region, the Ginzburg criterion is 
violated and mean-field theory breaks down. As the critical temperature is approached, 
low-energy fluctuations of the order parameter diverge owing to the flatness of the potential. 
Consequently, naive perturbation theory for loop corrections fails. One finds by dimensional 
analysis in terms of the correlation length £, that higher order interaction terms, e.g., the 
4-point coupling, diverge as the critical point is approached for d < 4, in particular in three 
dimensions. Therefore, a systematic analysis of the loop contributions in the critical region 
is in general difficult. 

In spite of these complications, various scaling relations have been found among the 
critical exponents. These relations imply that there are only a few independent critical 
exponents. In the O(N) theory, there are two independent exponents associated with the 
reduced temperature and the external field. Except for the hyper scaling relations, the 
scaling relations hold for empirically determined exponents in critical region, but also for the 
Landau mean-field theory. Scaling relations are easily derived, once the general assumption 
of homogeneity is made for the singular part of the thermodynamic potential density: 

F s (t,h) = L~ d F s (L A %L A *h), (8) 

where L is an arbitrary number not much greater than unity, d is the number of spatial 
dimensions, and A^h is the scaling dimension of the reduced temperature t and the external 
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field, h. This hypothesis was established more rigorously by Kadanoff using block spin 
transformations for the Ising model. Later on Wilson developed a systematic method, 
applicable to any system, for evaluating scaling dimensions A^h explicitly. The latter is 
known as the renormalization group (RG) method with the epsilon expansion about the 



critical dimension 



47 



48|. 



The RG method consists of two steps: i) integrating out a high momentum shell 
A/b < k < A with the parameter b > 1, and ii) rescaling the unit length k —> bk and 
other variables accordingly. Consecutive implementation of these procedures yields a flow of 
the renormalized Hamiltonian (thermodynamic potential), i.e., a flow under the RG trans- 
formation in the full parameter space V. 

A critical point of a continuous phase transition corresponds to a fixed point of the flow, 
where the length scale £ goes to infinity. Let V* be a fixed point, and v = V — V* the 
deviation from it. Then the thermodynamic potential density can be written as 

F(V)=F(V*;{v}). (9) 

Consider a system at a point in the parameter space, which is not a fixed point. After a 
single renormalization step, we obtain 

F(V*; {v r }, K» b~ d F(V*; {b A ^v r }, {b A ^v ir }), (10) 

where the deviations {v} can be classified into relevant {v r } and irrelevant parameters {vi r } 
according to their scaling dimension. By definition, the relevant (irrelevant) parameters 
have positive (negative) scaling dimension A vr > (A vir < 0), and ones with vanishing 
scaling dimension are called marginal parameters. The factor b~ d in front of F stems from 
the rescaling and reflects the dimensionality of F. Repeating this procedure n times results 
in the substitution b — > b n . 

Note that in general, a constant term appears in F after the renormalization procedures. 
This term, which breaks homogeneity, originates from integrating out the higher momentum 
shells. However, since this term is non-singular, it can be dropped. The remainder obeys the 
homogeneous relation: F s (V*;{v r },{v ir }) ~ b~ d F s (V*; {b A ^v r }, {b Avir v ir }). The homoge- 
neous scaling relation (JSJ) holds only close to the critical point, in the so called scaling region. 
Here the irrelevant variables are very small, and can be put to zero, since b nAv%r 1. The 
singular behavior near the critical point is controlled only by the relevant parameters, and 
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various scaling relations are obtained naturally, provided the system is sufficiently close to 
the critical point. The relevant variables {v r } again, can be identified with the temperature 
and the magnetic field, v T \ oc t and v r2 oc h. 



B. Static critical phenomena 

We first discuss renormalization of the effective potential Q = —(lnZ)/V per volume to 
define the static properties of the O(N) scalar field theory in the low energy limit. The 
partition function is defined by Z = e ~ H with the dimensionless reduced Hamiltonian 
7-L. The static renormalization group aims at tracing the evolution of the coefficients in the 
Hamiltonian (jlj), under the RG transformation. 

The theory is defined with a finite ultraviolet cutoff A. This means that the O(N) scalar 
field theory is an effective one, which can be applied only at scales below A. Since the soft 
modes are treated explicitly, the theory possesses the correct infrared behavior. We follow 



48 



This involves 



the renormalization group procedure developed by Wilson and Kogut 
the two steps mentioned above: integration over the momentum shell A/b < k < A in loops 
corrections with a parameter b > 1, and rescaling the variables and fields 

x -> x/b, (11) 
A bA, (12) 
& -> (13) 

The scaling dimension of the order parameter field, a^, is determined by the requirement that 
the rescaling leaves the auto-correlation function of 0; unchanged, i.e., keeping the derivative 
term of <pi to be marginal: = ^(d — 2 + r)'). Here r/ is the anomalous dimension, not be 
confused with the shear viscosity, rj. 

This procedure provides an evolution of the system under successive changes of the length 
scale and decimation of shorter wavelength modes. This process generates all couplings 
including higher order ones allowed by the symmetry of the system. The theory approaches 
a low-energy effective theory for the long wavelength modes. 

After repeating the renormalization procedure / times, one obtains the well-known recur- 
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sion relations for the coefficients, to leading order in the coupling u, 

r l+1 = b d ~ 2a - [n + 2(N + 2)fi 4 «i {A 2 (1 - b~ 2 ) - 2n In b}] (14) 
Ul+1 = b d ~ 4 ^ui [1 -4(JV + 8)fi 4 Ujln&] . (15) 

The above relations are obtained for d = A — e dimensions. The factor = 2 1 ~ d 7r~ d / 2 T(d/2) 
originates from the solid angle integration in d dimensions, divided by (2n) d , with T(x) 
being the Gamma function. In the right hand side of the recursion relations, the factors of 
b with exponents stem from the rescaling, while the terms proportional to Q 4 arise in the 
decimation of shorter wavelength modes. These two contributions play a competitive role in 
the RG evolution. This makes an appearance of non-trivial fixed points possible. A simple 
dimensional analysis shows that interaction terms higher than quartic are irrelevant under 
the renormalization. The recursion relations in fact admit a non-trivial critical fixed point, 

" '= 40^+8) +0 ' £j »' (17) 
implying that the system undergoes a second-order phase transition with infinite correlation 
length £. 

One can extract the scaling dimensions by observing how the coupling parameters behave 
near the fixed point. To do this, it is sufficient to linearize the recursion relations in terms 
of Sri = (r* — n) /f2 4 A 2 and 8u\ = u* — uf 



Sr l+1 -Sn \ / 2-f±| £ 4(A + 2) L , 

5u i+ i -5ui ) \ -e 




(18) 



Then eigenvalue problem of the above matrix tells that only r is the relevant parameter and 
Sr oc b Ar with A r = 2 — ^r§£ + 0(e 2 ) being the scaling dimension, while u is irrelevant, 
with a negative scaling dimension A u = —e + 0(e 2 ). 

We also note that in the long- wave length limit, the self-interaction of the field 0j vanishes 
in d = 4 because u* ~ 0(e). Thus, the perturbative expansion in the coupling constant u is 
equivalent to an expansion in e. This expansion is valid near the fixed point in a dimension 
slightly below four. In dimensions higher than four, the fluctuation contribution to the 
renormalization of 4-point coupling is negligible, i.e. the mean-field description remains 
valid. The physical correspondence of dimensionful quantity r — r* with thermodynamic 
variables is introduced by hand, e.g., r — r* oc T — T c near the critical point. 

11 



Let us now examine the interaction term ^^cfjE. Since only this term provides the static 
coupling between order parameter and energy density fluctuations, its critical behavior is 
crucial in the subsequent analyses of critical dynamics. A system with a non-conserved order 
parameter coupled to the conserved energy was classified by Hohenberg and Halperin (see 



45[ and references therein), as model C. The recursion relations for 7 and Cq are given 



Ref. 

by 

Cf + \ = b d ~ 2aE C^ 1 [1 - 2Nv x In b] , (19) 

v l+l = b d -^vi [1 - 8(N + 2)n 4Ul In b - 2Nv t In b] , (20) 

where v = VL^qCq is the dimensionless three-point coupling, and d— 2a^ — a/v — a6(a)/u 
with a being the exponent of specific heat C ~ (T — T c )~ a ~ t~ a . Here 9{x) is the unit step 
function and a is the critical exponent of the specific heat. The fixed point of the coupling v 
is given by v* = a + 0(e 2 ), which vanishes for negative a. The sign and value of the critical 
exponent a depends on N and d, as noted above. We return to this point in the subsequent 
section. 

For a > 0, fluctuations of the energy can become critical, i.e. the corresponding mass (the 
inverse of the specific heat) vanishes at the critical point. Thus, also the critical dynamics 
may be affected by energy fluctuations. On the other hand, for a < the mass term remains 
finite and fluctuations of the energy do not affect the static critical properties of other 
variables. Nevertheless, since the order parameter always exhibits critical fluctuations at a 
second-order transition, it is possible that these fluctuations affect other variables through 
dynamical effects, like mode-mode couplings. 

IV. CRITICAL DYNAMICS 

A. The stochastic equation of motion 



To address the critical dynamics of a system, one needs the equations of motion [49h51|. 
The low energy and long-wave length dynamics in the critical region is dominated by slow 
modes, i.e. fluctuations of the order parameter and the conserved quantities. We describe 
such modes by fields Ai(t, x) varying in space and time, and introduce a stochastic equation 
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of motion to describe the dynamics of the fields. In the mixed Fourier representation 

d t Mt, k) = L lm (k) - [A h A m ] PB + 9 ( (t, k). (21) 

oA m (t, k) dA m {t,k) 

Here H = T-L({Ai}) is a reduced effective Hamiltonian, which is a functional of the slow 
modes, and e~ H is proportional to the probability for a particular configuration of the 
fields A\. The first term on the right side involves transport coefficients Li m (k), which are 
responsible for the damping of fluctuations. Hence, this term describes irreversible processes. 
Owing to this term and the noise term G/, the system eventually reaches an equilibrium state 
where 8%/ '5 A\ = 0. The noise term satisfies the fluctuation-dissipation relation, 

(6,(t, k)Q m (t', k')) = 2L lm (k)5(t - t')5(k - k'), (22) 

which is valid for Gaussian noise. The cross terms with / ^ m originate from a possible 
bilinear mixing among the variables, ~ AiA m , in the effective Hamiltonian %. 

The second term, with the Poisson bracket, [• • -} PB , yields non-linear interactions, the 



mode-mode couplings [52 



531 ] . These describe the non-dissipative (reversible) processes, 
which are responsible for the large amplitude collective fluctuations induced by the critical 
behavior of the order parameter. Consequently, this term contributes to the singularities, 
which define the critical dynamics. The mode-mode couplings are formulated in terms of 
the generators of the relevant symmetries, and thus preserve the invariances of the original 
equations of motions. 

The equation of motion can be derived from the Liouville equation by the projection 
method in the Markovian approximation under some reasonable assumptions. The deriva- 



tion is reviewed in appendix A. Further details can be found in Ref. 54j. The presence of 



the Poisson bracket implies that the equations were derived from Hamilton's equations of 



the classical theory, which is valid for slow modes (see also discussion in Ref. |55|). 
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B. Response functions and transport coefficients 



A set of stochastic equations of motion for the slow modes Ai = E, J, Q^} is obtained 
from Eq. f l2~Tj) . given the effective Hamiltonian constructed above 01]), 



dE ^ -, 2 5U 5H n 

sT = r ° v ^-* V£ '^ + 9 ' 



Qjk\ PB X/O + 
Vj'fc 



a/ 



5U 



dt 



(23) 
(24) 

(25) 



m 

:lpB 5Q 



90 [Qij, 



Ik 



ij i Yk\ p p 



-go^Qij ■ — - + 9 

5J 



(26) 



where go and go are the mode- mode couplings associated with the translation and O(N) 
symmetries and T = 1 — 5^ is the projection operator on the transverse direction. Fluc- 
tuations of the transverse momentum describe diffusive modes, while fluctuations of the 
longitudinal momentum coupled with energy fluctuation describe sound waves, which have 
a linear dispersion relation. The latter are not taken into account because the sound mode 
corresponds to fast dynamics, which does not affect late time evolution. The longitudinal 
momentum, however, has to be considered in order to address the critical behavior of the 
bulk viscosity. This will be discussed in section IV. 

The Poisson bracket between and 0, are deduced from the quantum commutation 
relations 



WhQjklpB — —^j^ki + ^k^ij, 

[Qij, Qki] PB = —Qji^ik + Qjk^u + QuSjk — QikSji- 

In line with ( )22l) . the noise term correlation functions satisfy 

(6i(t, x)6 j (t', x')) = 2X 5 ij 5 d (x - x')5(t - t'), 

(9 E (t, x)9 E (t', x')) = -2T V 2 5 d (x - x')5(t - t'), 

(9j(t,x)9j(t',x')) = -2 Vo lV 2 S d (x - x')S(t - t'), 

(%(*, x)9 kl (t\ x')) = -2II (StkSfl - 6 a S jk ) V 2 S d (x - x')5(t - t'). 



(27) 
(28) 



(29) 
(30) 
(31) 
(32) 
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The transport coefficients are obtained from low energy limit of the dynamical response 
functions. In frequency and momentum space, they are given by 

i dx;\k) 

A = 'fcu^^' (33) 
I = ili^ dX i {k \ (34) 

l - = i^P dX t {k \ (35) 

Tj k^O did 

where we use the short hand notation k = {a;, A;}. The limit is taken first with respect to 
frequency and then to momentum. The response functions are obtained from the solution 
of the stochastic equations of motion, after averaging over the noise (■ • • ) & , 

- (37) 

= (38) 



5Hj(k) 



XQii(k) = (SikSji - 5 u 5 jk ) ( 5 ® %3 ^} ) 



The source terms are put to zero, after the variations in the above equations. 



(40) 



C. Dynamical renormalization group and dynamic scaling 

In this section we investigate the fixed points of the stochastic equations of motion for 
the O(N) scalar field theory, using the dynamical renormalization group (DRG) **. We thus 
determine the universal properties of the critical dynamics. The procedure of the DRG is 
very similar to the static one: i) in loop corrections, the momentum shell A/6 < \k\ < A 
is integrated out, while the frequency is integrated over the whole domain of definition 
— oo < ko < oo, and ii) rescaling of all variables. A difference from the static renormalization 
is that there appears a frequency scale, and its scaling in length units is assumed to be 

** In analogy to the static case, a fixed point of the equations of motion is necessary to be able to sort out 
the critical dynamics. 
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uj — > b z u> where z is the dynamical critical exponent for the slowest mode. The rescaling 
factors in the O(N) scalar field theory follow from dimensional analysis of the stochastic 
equations of motion and of the effective Hamiltonian, "H, which is a dimensionless quantity, 



X 


-> 


r 1 ^ 


(41) 


A 


-> 


6A, 


(42) 


t 


->■ 


b~*t, 


(43) 




— >• 


b a *<p, 


(44) 


E 


-> 


b aE E, 


(45) 


J 


-> 


b aj J, 


(46) 


Q 


-> 


b a ®Q. 


(47) 



The exponents for the fields are cle = (d — a/u)/2, aj = d/2, clq = d/2, and as in the static 
case a<£ = (d — 2 + r/)/2. We have set the scaling dimension of Xo to zero. 

In evaluating the dynamical response function, we employ a loop expansion in terms of 
the deviation from the upper critical dimension e = 4 — d in the same way as in the static 
case. For instance, the response function of the order parameter is expressed as 



[A + (48) 



where S represents the loop corrections integrated over the momentum shell and the bare 
propagator is given by 

G<t>{k) = Ty . (49) 

—iu + A ( r + k 2 J 

A renormalized relaxation rate for the order parameter fluctuation, A, is derived from the 
response function \4> pHj) . This procedure corresponds to a single renormalization operation. 
Thus, the recursion relation for A reads 



l + ^(A,,r { ,.-- ;6) , (50) 



where S is a dimensionless loop function, and the overall rescaling factor in b can be deter- 
mined from the rescaling factors of the other variables using the equation of motion. The 
recursion relations for the remaining transport coefficients follow the same procedure. See 
the following sections and Appendix B for details. 
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V. RESULTS FOR TRANSPORT COEFFICIENTS 



A. Flow equation 



In the dynamical renormalization procedure presented in the previous sections, we derived 
a set of recursion relations for transport coefficients to one loop order 



X l+1 = V-^'M 
Y l+1 = b z ~ 2 -^T l 
Vi+i = b z ~ 2 7]i 



1 - 



\i + ryq x lX i (A/ + Ui/xi) 



9i 



AY l C^{V l C; 1 + T ll ) 2 



A 2 

(1 - 6" a ) 



a 2 A 2 



9? 



(51) 

(52) 

(53) 
(54) 



The corresponding relations for the mode-mode couplings and static coefficients read 



9i+i 

9i+i 
n-\ 

li+i 
Xi+i 



b z - 2+e/2 ~gu 

b d-2a-a Eli [i _ 4 (at + 2) U 2 Q A In b - 2N 1 fC l ] 
Xh 



(55) 
(56) 
(57) 
(58) 
(59) 



where and oe were defined above. Note, that the mode-mode couplings g and g exhibit 
only trivial scaling without loop corrections. This follows from Ward identities for the 
higher order response functions. This can be also understood from Galilean invariance and 
invariance of the equations of motion under O(N) rotations. Using these recursion relations, 
we find the fixed points of the equations of motion, and extract dynamical critical exponents 
and scaling relations. 
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In the continuum limit, b — > 1, the recursion relations yield the flow equations, 



dT 
dr] 
dll 



A 

r 

v 
n 



Av* 

Z-2+T)'- - 

1 + Wi 

a 3 f 3 
z-2 + 



+ {N 



oo 2 



v A 1 + uj 3 



21 + w 4 2 J 



(60) 
(61) 
(62) 
(63) 



where d = J? with 6=1 + 5. We introduce effective vertices for the mode-mode couplings: 



<95 



11 

All 



^4, fz 



^ FT fi 4 A 2 , and f A ~- vILx 



dfx = fx 

dh = h 

dfs = h 

dh = h 



-2 + e-rj' + 



Av* 
1 + wi 



(N-l) 



h 



1 + u 2 24 



e-r]' + 



Av* 



- (iV- 1 



/a 



-2 + £ + 2iVt;* - 



3 /s 



4 1 + co 3 24 

3 / 4 1, _ J_ ' 
£ 2I + W4 2 /2 24 Jl 



l + w 2 2 /2 2 1 + o; 4 



/l 



(64) 
(65) 
(66) 
(67) 



and ratios of the transport coefficients, U\ 



vc- 



u 2 



w 3 



and 6^4 = 3^ 



9o;i 
du 3 



LU 2 



UJ 3 



LO.l 



I OAT * , 3 ^ AV * 

-7] — 2JSv + -—. 1- 



4 1 + cu 3 1 + ujx 



(N 



h 



1 + UJ 2 



3 U 



2 1 + 2 

3 / 3 



/a 



Av* 



1 + wi 



(iV 



/ 2 



1 + w 2 



-2AV + 



2I + W4 2 J 



4 1 + w 3 24 
1 



/1 



24 



fx 



(68) 
(69) 
(70) 
(71) 



Here the static fixed point of the three point function v* = a/(2Nu) has been inserted. The 
three point vertex is of order e, with ajv = (A — N)e/(N + 8) + 0(e 2 ) near four dimensions. 

The flow equations for the mode-mode couplings show that, except for f 2 , there are 
contributions of order 0(e°) on the right hand side. Since these equations admit only trivial 
stable fixed points, i.e. /* = f 3 = fl — 0, these mode-mode couplings vanish in the long 
wavelength limit. 

In the classification of the dynamical universality class we must, as implied by the discus- 
sion above, consider two cases, depending on the sign of a. The critical number iV c , where 
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a changes sign, is given by Fischer j47|: a is positive for N < N c , with N c ~ 4(1 — e) near 
four dimensions, and N r ~ 1.8 for d — 3. 



B. Fixed point for N = 1 



We first consider the N = 1 case, where the symmetry is reduced to the discrete Z2 
symmetry, and the energy fluctuation must be taken into account, owing to the small but 
positive exponent a > 0. The fixed points can be found by setting the right hand side of 
the flow equations to zero. To leading order in e we find: 

_ 2a 1 



A 



v 



2-2 



V 1 + UJi 

0. 



0. 



which admit a stable fixed point, 



u{ = 1, 



(72) 
(73) 
(74) 



a 

2 + -. 
v 



(75) 
(76) 



The last equation defines the dynamical critical exponent, which was deduced from the 
condition that T and A each have a finite non-trivial fixed point. Thus, the long-wavelength 
dynamics of the system is, up to order e, governed by fluctuations of the energy and the 
order parameter on equal footing. In the critical limit £ — > 00 (keeping finite) we obtain 
the following relaxation rates 



8<j>{t) ~ expi-Xx^kyH) ~ exp(-A; 2+a ^), 
8E(t) ~ exp(-rC~ 1 A; 2 t) ~ exp(-A; 2+Q/ ^), 



(77) 
(78) 



where we have used the fact that A ~ £ 2 z v ' ~ ^ a l v +°( £2 ) ; the order parameter susceptibil- 
ity X<f>{k) ~ k~ 2+v> , and T ~ ^ ^o+o(e 2 )_ result^ can also be obtained from the 



tt The dependence of transport coefficients (and of other physical quantities) on the coherence length (tem- 
perature) are derived as follows: let {a} be the full set of parameters (static coefficients) including relevant 
and irrelevant ones. Now we pick up only one relevant parameter a\, e.g., the reduced temperature a\ oc t, 
and set the other relevant parameters on the critical surface, i.e., to zero. Since £ = £({a}), a transport 
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fixed point of co^, which is the ratio of these two fluctuating modes, i.e., u\ ~ YC" 1 /X ~ £ v . 
Thus, to leading order in e, the critical exponent of A is smaller than that of T by a/ v. 

It follows from the discussion above that, owing to the dominance of the fluctuations of 
the non-conserved order parameter and the conserved energy, the single component scalar 



theory belongs to the dynamic universality class of model C. In Ref. [55[, the same conclusion 
was drawn based on the solution of a classical relativistic 4 theory on the lattice in d = 2 
spatial dimensions. 

An important point, which was not discussed so far is the renormalization flow of the 
shear viscosity. In the regime where z = 2 + a/u, rj does not reach a finite stable fixed point. 
This means that the critical dynamics of the order-parameter does not affect the shear 
fluctuations. Only short wavelength processes (rapid processes) contribute. Consequently, 
the shear viscosity remains finite, in contrast to model H, where a finite fixed point of a mode- 
mode coupling provides the scaling relation between the exponents of heat conductivity and 
shear viscosity. 

In order to obtain a finite fixed point of the flow equation for r]*, drj — rj [z — 2], we have 
to set 2 = 2, which is smaller than that of the fluctuating modes of the order parameter 
and the energy . Thus, long-wavelength fluctuation of the transverse momentum diffuses 
faster than the other modes, since 8J(t) ~ exp(—r] x^ 1 k 2 t) ~ exp(— k 2 t), where we have 
used a bare shear viscosity rj ~ £° and susceptibility Xri ~ £°- Therefore, fluctuations of 
the transverse momentum correspond to a faster mode and decouple in the long- wavelength 
dynamics inside the critical region. 

C. Fixed point for N > 1 

For N > 1 the static coupling between the energy density and the order parameter 



vanishes at the fixed point, as shown by Hohenberg and Halperin 45| (more precisely the 



coefficient T = T({a}) = T (£, {a}), where {a} represents the irrelevant parameters. We drop the ir- 
relevant parameters assuming that the system is very close to the critical point, and that the RG flow 
is sufficiently developed so the the parameters are in the immediate vicinity of the corresponding fixed 
point. Then, an RG transformation changes £ — > £/b and T — >• b x T. One thus finds the scaling relation 
b x T({a}) = T U/b, {b A a}) which leads to T - £~ x . 
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effective three-body coupling 7oC*o vanishes)". Therefore, the critical fluctuations of the 
order parameter do not directly affect the energy-momentum dynamics in the static case. In 
the dynamic case, such a coupling could be induced by the mode-mode coupling fi, which, 
however, vanishes in the long- wavelength limit. Moreover, critical fluctuations of the order 
parameter couple to the O(N) charge density only via the mode-mode coupling f 2 . Thus, 
for N > 1 one expects the energy modes to be irrelevant, while the O(N) charge fluctuations 
affect the critical dynamics owing to the mode-mode coupling f 2 . Taking these arguments 
into account, we find the fixed points of the flow equations in the same way as for N = 1 
case: 

N- 1 



A 

n 

h 

^2/2 



2 - 
- 2 



1 + oj 2 
1 



h 



0. 



0. 



N 



e — 



— f'2 - 77/2 
1 + co 2 2 



0. 



N-l 



0. 



i + w 2 : 

These equations yield the following stable fixed point and dynamical exponent, 



Wn 



2(N 



d 



(79) 
(80) 
(81) 
(82) 

(83) 
(84) 
(85) 



2 2 

The dynamical exponent is obtained by requiring that A and II have a non-trivial fixed point. 
At the critical point, the transport coefficients scale as A ~ £ £//2 and II ~ £ £ / 2 to leading 
order in e. These results are consistent with the fixed point of the mode-mode coupling: 
f2 ~ i~ e+v ' ■ Long- wavelength fluctuations of the order parameter and of the O(N) charge 
fall off with a characteristic frequency Uk ~ k d l 2 . 

Fluctuations of energy and transverse momentum are governed by the flow equations 
drj = 7](z — 2) and dT = T(z — 2). The fixed point at z = 2 implies that these fluctuations 



** As explained earlier, the critical dynamics is governed by the sign of the critical exponent, a. The absolute 
value of a is small for not too large N. Consequently, the sign of a is very sensitive to the approximations 
used. It is well known, that, to leading order, the epsilon expansion results in spurious sign of a in the 
range 2 < N < 4 [56l |. To obtain a physically correct result, we use the input from non-pcrturbativc 



methods according to which a in d = 3 is positive only for N=T, and negative otherwise. 
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are slower than those of the order-parameter and the O(N) charge with z = d/2<2ind<A 
dimensions. However, from the fixed point analysis we see that the critical fluctuations of 
the order parameter do not affect the energy and transverse momentum fluctuations in 
the long wavelength limit. Thus, although they participate in the critical dynamics at finite 
wavelengths, they decouple at late times. Consequently, owing to the dominance of the O(N) 
charge fluctuations the critical dynamics of the multicomponent O(N) theory is described 
by the dynamical universality class of model G. 



D. Bulk viscosity 



Before summarizing the main result of this work we briefly discuss the behavior of the 
bulk viscosity at the phase transition. The properties of the bulk viscosity in a slowly 



relaxing 



fluid and its possible singular behavior were first addressed in Ref. [5^] (see also 



Ref. 58j). The behavior of the bulk viscosity in system with a single component non- 
conserved order parameter was considered in Ref. 59|. Here, however, the critical exponent 
was not evaluated, but rather it was guessed based on input from experiment. 

In contrast to the shear viscosity, the bulk viscosity can diverge at the critical point in the 
O(N) model depending on the value of N. For the case of the single component scalar theory 
in d = 4 — e spatial dimensions, the bulk viscosity tends to infinity as £ ~ ^z-a/u _ ^2 ^ Q 
leading order in e), while for the multicomponent N > 1 theory the bulk viscosity remains 
finite ( ~ £°. Here z is the dynamical critical exponent, which was determined from the 
slowest mode function of N. 

In order to address the critical behavior of the bulk viscosity, the longitudinal component 
of the momentum current has to be considered in Eq. ( 125]) . In this case the projection 
operator on the transverse direction is dropped and an additional contribution owing to 
the bulk viscosity is added on the right hand side of Eq. f[2"5"j) . The critical behavior of the 
bulk viscosity can be deduced along the lines discussed in Ref. 
critical exponent for the bulk viscosity in model H was computed. Also, the QCD critical 
end point, which is theoretically expected to exist at a finite density and temperature in the 
QCD phase diagram 6jJ, belongs to the universality class of model H 62J. Recently the 
critical dynamics of the QCD critical end point was examined in a comprehensive manner 
based on DRG [63 1. 



60], where the dynamical 
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The results of Ref . [60| for the bulk viscosity can be immediately generalized to the single 
component scalar field theory since the (non) conservation of the order parameter does not 
affect the result as soon as the dynamical critical exponent is defined. Consequently, for 
N = 1 the bulk viscosity diverges jat the critical point as £ ~ ^z-a/u _ j nc i eec i ; the bulk 
viscosity is given by (see e.g. Refs 
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C = J^lim J™ dt J d d xe-^ (II«(af, *)%((), 0)) , (86) 

where is the stress tensor, which can be defined by comparing Eq. (j25p with the Euler 
equation dJi/dt = VjIIy. We are interested in only the dominant singular contribution 



to Eq. (I86p . As noted in Ref. 60|, it arises from the part of the stress tensor that is 
proportional to 70 2 . The integral in Eq. ( 186]) is taken over the domain with characteristic 
spatial extension of order £ and in the time direction of order C, z . Therefore, the dominant 
singular contribution to the bulk viscosity reads £ ~ £, z ~ d 1 2 X% which reduces to £ ~ 
after substitution of the renormalized quantities for 7 and \<t>- m contrast to model H, one 
should, however, keep in mind that in this expression z = 2 + ajv. Thus, extrapolating to 
e — > 1, we find that the singularity of the bulk viscosity is given by ( ~ £ 2 , while in model 
H it is stronger, ( ~ £ 2 8 . Note, that in both cases the ratio of the singular part of the bulk 
viscosity to the relaxation time of the O(N) charge fluctuations, r, vanishes at the critical 



point as £/r ~ £- a / u C" 1 , in agreement with 14|. This is a consequence of the fact that 
the single component scalar field theory belongs to the same static universality class as the 
liquid-gas phase transition. 

For iV > 1 the above discussion does not apply because the energy fluctuation decouples 
from the order parameter in statics, i.e. 7* — > 0, as we found in Section III. Owing to the 
vanishing mode-mode coupling g* — > in the long wave limit, the critical fluctuations do 
not couple dynamically to the current Jj either. 

Therefore, the bulk viscosity is finite ( ~ £° at the critical point. In this case, the ratio 
of the bulk viscosity to the relaxation time vanishes as (/r ~ 0(£~ z ). 

VI. SUMMARY 

In this paper we have evaluated the critical exponents for the dynamics of the O(N) 
scalar field theory with all possible slow modes. We showed that for the case of the single 
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component theory its dynamical universality class reduces to model C. The dynamical critical 
exponent is given by z = 2 + ajv. On the other hand, for the multicomponent theory, 
the critical dynamics is dominated by O(N) charge fluctuations. This drives the critical 
exponent down to the value z — d/2 and the theory belongs to the dynamic universality 
class of model G . In both cases, N = 1 and N > 1, the shear viscosity remains finite at the 
critical point, while the bulk viscosity diverges for N — 1, and remains finite for N > 1. 

In QCD, the 0(4) chiral symmetry in the light quark sector is broken by the finite u and d 
quark masses. For high temperatures and small values of the chemical potential, the second- 
order phase transition is replaced by a crossover. Our results imply that the singular part of 
the shear and bulk viscosity remain finite also at the QCD phase transition. However, from 
the present analysis within the DRG, we cannot draw any conclusions on the behavior of the 
regular parts of the viscosities near a second-order or a crossover transition. This problem 



can on 
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y be addressed in more microscopic approaches based on QCD or QCD-like models 
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671 ] . or within the novel microscopic approach to critical dynamics, employing 



the conjectured gravity dual description of conformal field theories 681-170] . 
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Appendix A: Derivation of stochastic equation of motion 

In critical dynamics we are interested only in tracing the evolution of slow modes. From 
microscopic point of view even if we start with a set of exact equation of motions for the slow 
modes, they would be inevitably affected by all the other degrees of freedom (including fast 
modes). Slow modes also would mix after finite elapse time. Therefore, we need a method to 
extract time evolution of slow modes only from full microscopic equation of motion, where 
the other degrees of freedom are fairly incorporated. 



1. Master equation with projection 

We first derive a master equation, i.e., an equation of motion for the distribution function 
Qaif) — 5 (A(t) — a) = UiS (Ai(t) — ai) , which defines a probability distribution for the 
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macroscopic variable Ai(t) to take the value ai at the time moment t. We start with Liouville 
equation 

±A l (t) = iLA l (t), (Al) 

where since we are dealing with slow modes, the operator L is supposed to be Poisson bracket 
with the classical Hamiltonian, 

iLAi{t) = [H,Mt)] PB . (A2) 

We would like to split g a (t) into systematic and fluctuating parts. At initial time we 
start from a state defined by the slow variables. In general there is no a priori rule for the 
choice of slow variables. The integral of motion are, however, required to be included among 
slow modes. The slow variables at initial time Ai(0) will be rotated in Hilbert space by the 
Liouville operator exp(itL). This would take Ai(t) out of the subset of slow modes. By 
the systematic part of g a {t) we mean the amount of an overlap between initial and elapsed 
distributions at time t. Therefore, it is reasonable to define a projection onto initial state 
with equilibrium average (•■■}, 

Pg a (t) = J2(9a(t)g b (0))g b (0) (A3) 

b 

with general properties of the projection operator such as PP = P, P + P = 1 and 
(Pg a (t)Pgb(t)) = 0. The time evolution of g a (t) = Pg a {t) + Pga{t) is given by 

—g a (x,t) = -J2 — [v l (a)g a (x,t)} + J db J ds(iLF a (s); b)g b (x, t - s) + F a (t), (A4) 



where 



(iLAi (0)g a _,. r ,, n s > 

Wa(0)) 

F a (t) = ie upL PLg a (0)- (A5) 

In deriving the above equation, we have used the following decomposition of Liouville oper- 
ator: e ltL = e ttpL + i J* ds e isL PLe^ t ~ s ^ L , which can be verified by taking time derivative 
from both sides. 
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2. Prom Master equation to Langevin equation 



One can derive non-linear Langevin equation for Ai(t) from Eq. (1A4|) by taking the first 
moment of the distribution function, Ai(t) = J daaig a (t), 



d 

dl 



Ai(t) = v l [A{t)} + J ds(iLRi(8y,A(t-8)) + Ri(t), (A6) 
~ M,L(0)W;^)> + (W0)i ^ (0)}pfl . A(t)) ?^(*)] 



-^.[AWl^M + ^W, (A7) 

where we introduced the effective Hamiltonian for macroscopic variables H(a) = — ln(g a ) in 
units of /cgT = 1, and 

vi[A(t)] = ^(«)U=aw = (iLA(0);a)\ a=A{t) , (A8) 
i2j(t) = e itpL PM0) = ie itpL PLAi(0), (A9) 



oc 



L lm [A(t)] = y ds Uj~W) U=A{t). (AlO) 

Note that in derivation of Eq. (1A7|) from Eq. flA6[) the Markovian approximation for the 
memory term £ ds ■ ■ ■ A(t - s) -> J °° ds • • • A(t) was applied. We also assumed that the 
background transport coefficient Li m (A) is approximately independent on Ai(t) at late times. 
Owing to properties of Poisson brackets the first term in Eq. (1A7I) vanishes in most cases. 

Important point here is that since Ri(t) oc P, the force Ri(t) is uncorrelated with any 
macroscopic variables by construction (G[A(0)]Ri(t)) = for any arbitrary function GL4(0)]. 
In this sense the force Ri(t) is a pure random force. 

The second term in Eq. (1A7I) . known also as mode- mode coupling, describes reversible 
process, and involve non-linear interactions among Ai(t), responsible for critical dynamics. 
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3. Application to O(N) model 



Substituting the slow mode candidates in O(N) model and their effective Hamiltonian to 
the above Langevin equation, we obtain 



dt 

dE 
~dt 
df 
~dt 



-A (r (j)i - V 2 4>i + u o 0f 4>% + 2 7o0i-E - hA - g 
+2<?o0j (XqQh - /%) + 0i> 

T V 2 ((J^E - V 2 E + To 2 + p)-g (v#) ■ ( J - # 



J - H 



T 



dQ 



AB 



dt 



(All) 
(A12) 



r?oV 2 yJ - Hj + #o (V& J (r & - V 2 & + u o 2 0i + 2 7o <^£ - V 
+S V£ (c^E - V 2 E + 7o 2 + /?) + fjo^ Q ab {Xq'Qab ~ Hab) + 0j] , (A13) 
n V 2 {xq^ab ~ Vab) ~ 9o (^Qab) -(j-H) 

+9o U B h A ~ <pAh B ~ <Pa^ 2 <Pb + </>bV 2 (/> a ) + 6 AB . (A14) 
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In Fourier space k = {u, k}, the formal solution is given by 



AoMfc) + Oi(k) - ig / Hh) h -\j(k- k x ) - H{k - k x ) 



-Xqu / (j)(h) j (j)(k2)j^i(k -h- k 2 ) - 2A 7o / <t>i{k x )E{k - fci) 



12 



+2^0 / <Pj(h) [x Q 1 Qij(k - h) - Hij(k - kx)] 



(A15) 



E(k) = G° E (k) x 



T k 2 p(k) + - r o7 o^ 2 y <t>i{ki)<t>i{k - h) 

-ig j ' Eik^kx^Jik-k^-Hik-ky)} , (A16) 
= G° (*) • % ■ [ Vo k 2 H(k) + ((k) 

+igo J^ki<j)i{ki) jx^(fc- ki)(f>i{k - fci) - - fci) + 270^ <fii(k 2 )E(k - fci - k 2 ) 
+ig u / ki<f) i (k 1 )(f) j (k2)(l)j(k 3 )(l) i (k -k 1 -k 2 - k 3 ) 

J 123 

+ igoj^hEih) ^x E \k-ki) E ( k ~ k i) + P( k ~ k i) + 7o jj>i%)<l>i{k - k x - fc 2 )J 
+igo J QAB(h)ki [xQ~QAB(k - fci) - /iAe(fc - fci)] , (A17) 



Q AB (k) = G° Q (k) 



Il k 2 fi AB (k) + 6 Q {k) - ig / Q A B(ki)k! -{J(k- fci) -H{k-k 



+9o / {<pB{ki)h A (k - k x ) - (f) A (k 1 )h B (k - k x ) 
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(A18) 



where X,* (fc) — r o + fc 2 the inverse of static susceptibility, / 12 ... m — / n™ =1 du; n d d fc n /(2 
and (7fe)j- = 5jj — kikj/k 2 . Bare propagators read 



,d+l 



Cj(fc) ■ % 
G° E (k) 

Gl(k) 

G° Q (k) 



—ibj + r\k 2 



-Ik: 



-iu + r fc 2 ^Co" 1 + fc 2 
i 

-iuj + A fr + k 2 
1 



Note that G°Ak) and Gn(fc) are of diagonal form in O(N) space. 
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The noise-noise correlation functions satisfy the fluctuation-dispersion relations, 



0W?i(*O> 



mm) 



2A S i:i S [k + k'j 
2T k 2 s(k + k') , 
2r) k 2 S i:j S (k + 9 



) 



(A23) 
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The renormalized transport coefficients are determined from the response functions, 



Appendix B: Response function 

The response functions of slow modes are obtained from a set of stochastic equation 
of motions. The old-fashioned perturbation method, i.e., iteration of formal solution for 
different orders of interaction terms and taking average over noises in energy-momentum 
space, systematically generates loop corrections to response functions. In this article we 
perform only leading order calculations. In order to proceed with the calculations beyond 
the leading order, it is preferable to use the alternative field-theoretical approach (already 
mentioned in the text). 

1. Order parameter relaxation constant 

The response function for order parameter fluctuations is given as follows: 




(A26) 




(A27) 




(A28) 




(Bl) 



where loop corrections are accounted for in H^(k). 
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o 

V / 



(c) 



(d) 




FIG. 1: Leading order contributions to x<t>- Solid line represents propagator of order parameter fluc- 
tuations, dashed - energy, wavy - transverse momentum, and double - O(N) charge fluctuations, 
respectively. The solid circle indicates external field. 

The diagrams for the leading order contributions are shown in Fig. [TJ and the correspond- 
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ing equations are given by 



EJ(fc) = 2(N + 2)(\uY 
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1.2 



r + ( k — k\ — k2 
1 



r + k\ 



r + kl 



-iu + X 



3r + Ik - ki - k 2 ) + k\ + k 
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The renormalized order parameter relaxation constant to leading order is given by 



rcn 9 (-iw) 



A" 1 



fe->0 



4 7 2 CA 
A + T/C 



SI4 In 6 



ff 2 (iv-i) 

Axq (A + U/xq) 
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(B6) 



where fc± = k% ± We evaluate the above equation near d = A and at the critical point, 
where a renormalized mass goes like r ~ (note that non-trivial fixed point r* can be 
adjusted to 0). Since ~ k 2 , X acquires no mode- mode coupling contribution of order of 
g 2 at zero momentum. 

In the following, we take the same procedure to obtain the other transport coefficients 
renormalized to the leading order. 
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2. Energy diffusion constant 
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FIG. 2: Leading order loop corrections to xe- 



The response function for energy fluctuation is given by 



XE(k) 



I 5E(k) 



G E (k) Tk 2 + E E (k) 



Leading order contributions T, E are depicted in Fig. 2, and are given by 
E E (k) = G E {k) (2A 7 rP) 2 f G^kx) [G^-k^G^k - h) + G^k - fci)G^(fci 

2 r 2(r + k\ + P/A) l 



G E (k)N2\( 1 Tk _ 
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C - Tk 2 G E (k) 
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k-T % _-k 



-iu + TC~ l k\ + rjk 2 



where k± — k\ ± k/2. 

Evaluating equations above near d = 4 and at the critical point: 
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3. Shear viscosity 
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FIG. 3: Loop corrections to xj- 



The response function for the transverse momentum is of a tensor form because of the 
projection operator T, with leading order corrections £j shown in Fig. 3, 

5Ji(k) 



Xj(k) (TiV. 



5Hj(k) 



3 3 



(Bll) 
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;bi3) 



where k± — > — yk\ — k/2j to reach the last equality in £}. The transverse part Sj is 
extracted by projection (Tr7g = d — 1), 
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For r ~ 
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where (■ ■ - )a implies taking an average on the solid angle at d = 4: a d- dimensional angle 
average of an even power of one momentum component is defined by 



<*r>d = 



Jd d kk 2 -6[k 2 -l) !. 3 ...(2n-l) (d/2-1)! 
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For the present use {kf)± = \ and (kf)^ = |. The renormalized 77 eventually reads 
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4. O(N) charges diffusion constant 

The response function for the O(N) charges diffusion constant is defined by 
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FIG. 4: Loop contributions to xq- 
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where different contributions depicted in Fig. H] read 
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The renormalized O(N) charge diffusion constant is given by 
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where we have taken the solid angle average: 
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